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We present a comprehensive analysis of the relaxation dynamics of a Luttinger liquid subject to 
a sequence of sudden interaction quenches. We express the critical exponent f3 governing the decay 
of the steady-state propagator as an explicit functional of the switching protocol. At long distances 
ft depends only on the initial state while at short distances it is also history dependent. Continuous 
protocols of arbitrary complexity can be realized with infinitely long sequences. For quenches of 
finite duration we prove that there exist no protocol to bring the initial non-interacting system in the 
ground state of the Luttinger liquid. Nevertheless memory effects are washed out at short-distances. 
The adiabatic theorem is then investigated with ramp-switchings of increasing duration, and several 
analytic results for both the propagator and the excitation energy are derived. 



I. INTRODUCTION 

The on-going experimental activities on ultracold 
atoms are continuously challenging our understanding 
of many-body quantum systems. 1 Laser- or evaporative- 
cooled below the /iK, atoms crystallize in artificial 
lattices 2 thus providing nearly ideal realizations of 
bosonic 3-6 and fermionic 7-9 model Hamiltonians. The 
possibility of tuning the model parameters in real time 
brought the attention back to a fundamental issue in 
quantum statistical physics: does a non-interacting bulk 
system relax toward the correlated ground-state upon the 
switch-on of the interaction? 

"Sudden quench" is the nomenclature coined for the 
sudden change of a parameter like, e.g., the convexity of 
a parabolic trap or the interaction strength, in an equilib- 
rium system. ' n During the last five years experimental 
and theoretical investigations on the relaxation proper- 
ties of quenched ultracold atoms enlightened the intrigu- 
ing phenomenon of the thermalization breakdown: 12 ei- 
ther the system does not reach a steady state or, if it does, 
the steady state is not the ground state of the quenched 
Hamiltonian. A sufficient criterion for the occurrence of 
a steady-state has been found, so far, only for integrable 
models 13-16 and it has been argued that steady-state val- 
ues are calculable by averaging over a generalized, initial- 
state dependent Gibbs ensemble. 17 

The thermalization breakdown poses questions which 
are certainly conceptual in nature but may also be rele- 
vant to the growing field of optimal control theory: 18-20 
what is the steady-state dependence on the initial state? 
and on the switching protocol? how the adiabatic limit 
is recovered? In this paper we provide a comprehensive 
analysis of the behavior of a Luttinger Liquid (LL) sub- 
ject to arbitrary interaction quenches. We extend the 
study of Cazalilla for a sudden quench 21 to a sequence of 
N sudden quenches using a recently proposed recursive 
method. 22 Continuous quenches of duration T are then 
obtained in the limit W^oo and allow us to address the 



adiabatic limit by making T larger and larger. We cal- 
culate the equal-time one-particle propagator G^(x,t) 
as well as the excitation energy. At long and short dis- 
tances the steady-state G^ (x, t) ~ x~^ and in both 
cases we are able to write the critical exponent (3 as an 
explicit functional of the switching protocol. In the limit 
of continuous quenches (A — > oo) the propagator at short 
distances thermalizcs whereas at long distances does not. 
An analytic formula for ramp-like switchings of duration 
T valid for all x is derived and is shown that G^ and 
the ground state propagator are the same up to a critical 
distance that diverges for T — > oo. The recovery of the 
adiabatic limit is further illustrated from energy balance 
considerations. The calculation of the excitation energy 
AE(T) is reduced to the solution of a simple differential 
equation that can be used to find the optimal switch- 
ing protocol of duration T that minimizes AE{T). We 
prove that AE(T) > (strictly positive) for all switch- 
ing protocols of finite duration and provide an analytic 
expression for ramp-like switchings. 

The paper is organized as follows. In the next Section 
we introduce the model and the recursive procedure to 
calculate G and AE. Results for arbitrary sequences 
of sudden quenches are here illustrated. The limit of 
continuous switching protocol is carried on in Section III 
along with the derivation of several analytic formulas. A 
summary of the main findings is finally drawn in Section 
IV. 



II. SEQUENTIAL QUENCH IN A LUTTINGER 
LIQUID 

The sudden interaction quench in a LL has been ad- 
dressed in a series of papers. 21,23-28 At the distance x 
the propagator G^(x,t) exhibits the "light-cone" effect, 
i.e., a crossover between Fermi liquid behavior for times 
t <C x/2v (v being the quasiparticle velocity) and non- 
thermal LL behavior in the long time limit. Sequential 
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quenches yield an even richer phenomenology since addi- 
tional time (and hence length) scales appear in the prob- 
lem. We will show that different steady-state regimes 
emerge by probing the system at distances shorter or 
longer than the quenching time (in units of v) and that 
their nature depends on the switching protocol. 

The LL Hamiltonian describes interacting spinless elec- 
trons confined in a ID wire of length C and reads 

1 ^ f c/2 

+ 9iP 2 a (x) + 92P a (x)pa(x)], (1) 

where a denotes the chirality of the electrons with Fermi 
velocity e a v F {e R/L = ±1) and density p a =: ip^ipa :, 
": :" being the normal ordering. The coupling constants 
.92(4) refer to forward scattering processes between elec- 
trons of opposite(identical) chirality. We consider the 
system noninteracting (g2 — 54 = 0) and in the ground 
state before a series of N interaction quenches at times 
= to < t\, ■ ■ ■ < tjf-i = T takes place. Let j m be the 
value of the couplings gi—2,4 between t n -\ and t n and 
H n the corresponding LL Hamiltonian. Each H n can 
be bosonized 29 in terms of the scalar fields 4> a defined 
from ib a (x) = , K ° e »y^0 a (x) w ^ K the anticommut- 

v ' \2-na 

ing Klein factors and a a short-distance cutoff. The result 
is a simple quadratic form 

v f°° 

H n = ^- dx[K~\d x ct>{x)) 2 + K n (d x 6(x)f], (2) 

where <j) = 4>r + 4>l an d 9 = <pR — 4>l ar c con- 
jugated fields, v n = ^/{2nv F + g in ) 2 - gj„)/2n is 
the renormalizcd velocity and the parameter K n = 
yj (2ttvf + 9An - 92n)/ {2-kvf + g\n + 92n) measures the 
interaction strength. Note that < K n < 1 for repul- 
sive interactions; K n = 1 corresponds to noninteracting 
systems while small values of K n indicate a strongly cor- 
related regime. In our case Kq = 1 but no complications 
arise from arbitrary values of Kq, which is therefore left 
unspecified. As we shall see, such freedom permits to 
address general initial-state dependences. 

A. Excitation energy 

The excitation energy AE(t) is defined as the differ- 
ence between the energy of the LL at time t and the 
ground-state energy of the LL Hamiltonian at the same 
time. Since we are interested in AE(t) after the inter- 
action quench is completed we consider t > tj^-i and 
write 

AE(t) = (* (t)\H N \9 (t)) - (* N \H N \V N }, (3) 

where |* (*)) = e"^^"^- 1 ) . . . e -iHl (* I-t °)|tfo) and 
here and in the following \^f n ) is the ground-state of H n . 
To calculate AE(t) we expand the scalar fields in H n as 



where J\f a is the number of electrons with chirality a. 
Then, the Hamiltonian takes the diagonal form 

H n = J2Y1 Vnq b Un b °"in + E n , ( 4 ) 
q>0 a 

with E n the zero point energy and baqn the annihilation 
(creation) operators for elementary excitations of chiral- 
ity a and momentum q of H n . These operators are re- 
lated to the non-interacting baq of the cj) a expansion via 
the Bolgoliubov transformations b aqn = b aq cosh ip n + 
bl q smhip ni with^„ = itanh- 1 [(l-A'2)/(l + J fC2)]. Fol- 
lowing the bosonization the average over the ground state 
\^ n ) is converted into an average over the vacuum of the 
^^-excitations and hence AE(t) takes the form 

AE(t) =Y,vqn q (t), (5) 

where n q (t) is the average of the number operator 

J2 a b a q N-i b aqN-i over |*o(<))- The relation between 
two consecutive boson operators are easily found and 
read 

baq(n+l) = b a qn 

cosh Atp n + sinh A(p„, (6) 

where Atp n = ip n +i — Vn- After a cascade of transforma- 
tions (6) to express the b aqn in terms of b aq o we obtain 
n q (t) = p q fl as the solution of the recursive system of 
equations 

P g ,n = (<£ + sl)P q , n+1 + 2c n s n Re[Q^ n+ ie- 4£ ' At "], 

Qq,n 2c n S n P q ^ n -\-\ 

+ r 2 (l M „-!£,Ai„ , 2nt i£,At„ 

P q ,n = Pq.n+i + 2s^P g , n+ i + 2c„s n Re[Q, in+ ie~ ie?At ™] , 

(7) 

with At n = t n+ i - t n , t N = t, c n = coshA(p„, s n = 
sinhA(^„, e q = 2vq and boundary conditions P q ^ = 1, 
Qq,N = Pq.N = 0. In Eq. (7) we assumed that v n = v for 
all n; the general recursive scheme is simply obtained by 
replacing vt n — > v n t n . 

In Fig. 1 we plot n q (t) as a function of q (for t > 
tftr_i = T there is no dependence on t since the switching 
protocol finishes at T) for a series of N = 4, 6, 7 and 00 
quenches at times t n = nT/(N — 1) with K n = 1 — 
0.35n/A and T = 1. The distribution of g-excitations is 
peaked at q = but for any finite N there is a revival 
every time g/[27r(AT—l)] is an integer. In the next Section 
we derive an analytic result for N — > 00 and show that n q 
vanishes exponentially at large q. As n q is not identically 
zero the LL does not thermalize. Below we show how the 
thermalization breakdown is reflected on the equal-time 
one-particle propagator. 
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FIG. 1. Excitation number n q versus q for N = 4, 6, 7 and oo 
quenches at times t n = nT/(N — 1) with _R" n = 1 — 0.35n/iV 
and T = 1. Momentum q is in units of 1/a and time in units 
of a/w. 
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FIG. 2. Density plot of In \G [N] (x, t)\ (arbitrary units) for 
JV = 3 as a function of time and distance. The parameters 
are K = 1, Ki = 0.1, A' 2 = 0.5, K 3 = 0.9 and T = 0, Ti = 
5, T2 = 10. Time is in units of lOa/u and distance is in units 
of 10a. 



B. The equal-time propagator 

The equal-time one-particle propagator is defined as 

G^(x,t) = (*o(*)|^(a04(O)l*o(t)>, (8) 

where the superscript N specifies the number of sudden 
quenches of the switching protocol. To calculate G^ we 
express the fermion fields in terms of the boson fields, ex- 
pand the latter in elementary excitations and exploit the 
transformations (6) between two consecutive b operators. 
We then obtain 

J_ p T. q>0 ^ aq T- q {^^<i^\Fn q ,o\ 2 ~\FL q , Q \ 2 ) 



2ira 



GW(x,t) 

(9) 

where the functions F aQt o(x, t) are solutions of the recur- 
sive relations 



otg, 



q.n 



± aq.n+1 c <--ti 



ivqA.tn 



(10) 



with boundary conditions Fn q ,N = (e lqx — l)cosh<£>jv, 
Fh q ,N = {e~ lqx — l)sinh< ( 9Ar. Like in Eq. (7) the general 
recursive scheme is obtain by replacing vt n — > v n t n . 

The full analytic expression of G^ grows in complex- 
ity with increasing N. To illustrate the typical features of 
qW] in Fig. 2 we report the contour plot of In \G^ (x, t)\ 
for N — 3. The propagator has a Fermi liquid exponent 
in the light-cone region t < x/2v, in agreement with Refs. 
21 and 30. At the quenching times t n the propagator ex- 
hibits a time-derivative discontinuity due to the sudden 
change of the interaction which, being a global perturba- 
tion, instantaneously affects the whole system. The n-th 
quench gives rise to incoherent excitations that propa- 
gate at a speed v and those at distance x scatter at time 



t = t n + x/2v. When x = L mn = 2v(t m — t n ) the scat- 
tering time coincides with the time of the m-th quench 
and a pronounced peak in In \G^- N \x, f)| as a function of 
x develops and persists forever. The peaks are visible as 
horizontal lines in the contour plot at x = Liq = L21 and 
L-2{). Besides the primary light-cone patterns with ori- 
gin in x — and t = t n we can see secondary light-cone 
patterns with origin in x = L mn and t = t m which in 
turn generate ternary light-cone patterns and so on and 
so forth in a cascade that grows like 2 N . 

In Fig. 3 we display the steady-state value 
lim t _ i . 00 G^(x,t) with spikes at the characteristic length 
scales L mn , in accordance with our previous discussion. 
The analytic expression of limt_ i . 00 (x, t) is not a 
simple power-law. Nevertheless, a power-law behavior 
~ x _/3 is recovered at short and long distances (see inset 
of Fig. 3). By employing the recursive method of Eq. 
(10) we found for the steady-state behavior of G>- N ' at 
short distances 



lim lim G [N] (x,t) w — ; 

x-»0 t-s-oo 2tt[x + ia) 



i[K n 



(11) 



with a history and initial-state dependent (through Kq) 
exponent 



1 



K 2 

N 



2Kn 



N 1 



1 



K n 



-1. (12) 



The product structure of this result is similar to that 
of the exponent /3 in the TV characteristic / oc V 13 of an 
out-of-equilibrium LL subject to a sequence of interaction 
quenches. 22 Note that /3 s d only depends on the interac- 
tion parameters K n and not on the switching times t„. 31 
Equation (12) returns the well-known exponent of a LL in 
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FIG. 3. Steady-state value of |G' 3 ' versus x for Kq = 1, K\ = 
0.8, K 2 = 0.6, K 3 = 0.4 and T = 0, T x = 10, T 2 = 20. 
The crossover between the short and long distance regimes is 
clearly visible in the inset. Distance is in units of 10a and 
G [3] is in units of 10~ 4 /27ra. 



the ground state for Kq = K\ . . . = Kn = K since /3 s d = 
(K + 1/ K) /2 — 1 , and the exponent of a LL after a single 
interaction quench for Kq = K\ = . . . = K n = 1 and 
K n+1 = ... = K N = K since /? sd = (1 4- K 2 ) 2 /(4/i 2 ) - 1, 
in agreement with Ref. 21. 

The situation is radically different at long distances. In 
this limit the dependence on x of the steady-state prop- 
agator is again a power-law 



lim lim G [N] (x,t) 



>oo t— foo 



27r(x 4- ia) 



R[t n , K n ] 



MK n ] 



(13) 

but i?[t n ,/£T n ] is a length depending on all intermediate 
switching times t n and interactions K n . As for the ex- 
ponent /3m a striking cancellation of the intermediate K n 
occurs and we find 



Pu[K n 



K) 



AKqK 2 n 



(14) 



which depends only on the initial state. 

To summarize the steady-state propagator does not 
thermalize for discontinuous switchings, and the thermal- 
ization breakdown manifests in different ways at short 
and long distances. In the next Section we address the 
evolution of this behavior when the sequential protocol 
approaches a continuous protocol. 



III. CONTINUOUS QUENCHES 

Continuous switching protocols arc obtained as a lim- 
iting case of a sequential quench. Let us start by an- 
alyzing again the excitation energy AE(t). Taking the 



quenching times t n = nt/(N — 1) equally spaced and 
letting N — > oo the variable t n becomes a continuous 
variable between and t and we can construct the dif- 
ferentiable functions P q (s), Q q (s) and tp(s) according to 
P q {t n ) = P q ,n, Qq(t n ) = Qq,n and <p(t„) = ip n . Expand- 
ing Eqs. (7) to first order in At n = At = t/(N - 1) and 
Aip n we find a coupled system of differential equations 



dP q 
ds 

dQ q 
ds 



-2MQq]% 



-2P^ 
ZFq ds 



2ivqQ 



(15) 



that should be solved with boundary conditions P q {t) = 1 
and Q q {t) = 0. The average occupation number n q (t) 
is then simply given by 1 — P q (0). The most popular 
continuous protocol is the ramp-switching 33 ip(s) = (ps/t 
with < s < t. In this case the system (15) can be solved 
exactly and we find 



-*>- ^-T-V'" (16) 

which correctly approaches zero for t — > oo (adiabatic 
limit), whereas for any finite t is exponentially suppressed 
at large q, see the curve N — > oo in Fig. 1. We also 
checked that inserting Eq. (16) into Eq. (5) and taking 
the long-time limit t — > oo and the weak-coupling limit 
tp <C va/t the excitation energy vanishes as 



. „ , x „ In t 
lim AE(t) ~ L^r, 



tp <C va/t, 



(17) 



in agreement with Ref. 28. It is interesting to observe 
that the dependence of AE ~ 5 2 In 5 on the ramp rate 
5 = (p/t is not a simple power- law 8 V and, therefore, does 
not belong to the non-analytic regimes contemplated in 
Ref. 34. 

The system of differential equations (15) is exact (non- 
perturbative) and we now exploit it to prove that there 
exist no switching protocol of finite duration capable to 
drive the initial non-interacting system in the interacting 
ground state of a LL. Since AE(t) is the sum of non- 
negative n 9 's it is sufficient to show that n (t) = 1 — P (0) 
cannot be zero. The optimal switching protocol tp(s) that 
reproduces a target Pq{s) follows from Eqs. (15) with 
q = 0; in this case <p(s) depends only on the instantaneous 
value of Pq(s), and reads 



1 



ip(s) = <p — — cosh Po(s), 



(18) 



where (p is the value of the correlation angle at the end 
of the quench. Thus, -Po(O) = 1 only provided that the 
initial and final strength of the interaction is the same, 
<p(0) = (p. In particular no (t) > for all switching pro- 
tocols that connect an initial interaction strength Kq to 
a final interaction strength K ^- Kq. 

Analytic results can be obtained for the one-particle 
propagator as well. Let us start by analyzing the ex- 
ponents j3 of the long- and short-distance power-law 
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FIG. 4. Steady-state propagator |G'°°'| versus x for two dif- 
ferent ramp switching-times (in units of 10 2 a/v) T = 0.1 (blue 
solid) and T = 20 (red solid) and tp — 0.46. For comparison 
the equilibrium propagator with exponent fi s d (dotted) and 
the sudden quench propagator with exponent /3\d (dashed) 
are also shown. Distances are in units of 10 2 a and |G'°°'| is 
in units of 10 _2 /a. 

behavior of G^ N \ For a continuous switching we can 
construct the diffcrentiable function K(s) according to 
K(t n ) = K n , with K(0) = K and K{t) = K. Approx- 
imating K(t n -i) = K(t n - ^) w K(t n ) - 
taking the logarithm of Eq. (12) we find 



K'(t n ) and 



In 



2K (p s , 



1] 



lim 



N- 



N 

X 

n=l 



loe 



1 



T K'(t r 



and hence 



Psd[K(s)] 



K 



1 

K 



= ln 



1, 



N K{t n ) 
K 



K 



(19) 



which coincides with the exponent of a ground-state LL. 
Therefore, at short distances the initial-state dependence 
as well as the history dependence are washed out by dif- 
fcrentiable switchings and the propagator thermalizes. 
Discontinuous switchings 32 do, instead, introduce a his- 
tory dependence through the ratio between the values 
of K(s) across the discontinuity. To the contrary the 
long-distance exponent /?id in Eq. (14) is always history 
independent and the steady-state propagator docs never 
thermalize. Our conclusions agree with recent perturba- 
tive results by Dora et al. in Ref. 28. 

To study the crossover from short to long distances we 
must calculate the steady-state propagator at all x. For 
continuous switching the recursive system of Eq. (10) 
reduces to a system of differential equations 

dtp 



±F 

ds aq 



ivqF aq + F£ 



ds 



(20) 



that should be solved with boundary conditions Fn q (t) = 
(e iqx -l) cosh tp(t) and F Lq {t) = (e~^ x -l) sinh<p(i). The 
equal-time propagator G^ 00 ' (x, t) can then be calculated 
from Eq. (9) with F aq $ = F aq (Q). In the special case of 
a ramp protocol ip(s) = <ps/T for s < T and tp(s) = Cp 
for T < s < t the system of equations (20) can be solved 
analytically and we get 

(l- e «')sinh 7g [eivq(t _ T) _ ^ 
+e -^ t - T \t q vT- lq ^-±^)coshg,}, (21) 



where 7 9 = \/(p 2 — (vqT) 2 ; the function F£ (0) is ob- 
tained from Fnq(0) simply by exchanging sinht^ o 
coshi^ and by replacing v — > — v. In Fig. 4 we plot 
the steady-state propagator |G'°°1| as a function of the 
distance x for two different ramp switching-times and 
Cp = 0.46, which corresponds to the LL parameter K = 
0.4. The results clearly agree with the scenario outlined 
above: at distances smaller than the characteristic length 
2vT the steady-state propagator behaves thermally while 
at large distances it behaves like the sudden-quench (non- 
thermal) propagator of Ref. 21, albeit shifted upwards 
by a history-dependent constant (see also Ref. 28). 

The long-to-short distance crossover is particularly 
transparent in the weak coupling limit tp -C 1. In this 
case we can expand -Fq^O) to second order in Cp, perform 
the sum over q and find the steady-state propagator 
at all x 



lim G [ooI (a;,t) 



2-k(x + ia) 



x 2 - (2vT) 2 



2vT 



x 2 


x tp 
4v' 2 T 2 


x - 2vT 


x 2 - {2vT) 2 




x + 2vT 



(22) 



with power-law exponents depending on x. Equation (22) 
reproduces with remarkable accuracy the long-to-short 
distance crossover, which in the same approximation read 



lim hm GW(i,t) W . 

ir->0 t-voo 2ir(x + ia) 



lim lim G [oai (x,t) w — , — 

i^oof^oo 2tt{x + ia) 



29- 



(23) 
(24) 



R = \2vTa being a history-dependent length [cfr. 
Eq.(13)]. The high accuracy stems from the fact that 
the correlation angle (p remains small also in the inter- 
mediate coupling regime characterized by gi w 54 w Dp 
for which tp w 0.1. 

Finally we observe that from Fig. 4 and Eq. (22) 
the adiabatic limit is recovered when the switching time 
T — > 00. As T increases the steady-state propagator 
equals the equilibrium propagator up to larger and larger 
values of x. 



6 



IV. CONCLUSIONS 



We provided a comprehensive analysis of the relaxation 
dynamics of a LL after the quenching of the electron- 
electron interaction for different switching protocols. The 
bosonization method is combined with a recursive proce- 
dure to address arbitrary sequences of sudden quenches 
and hence, as a limiting case, continuous protocols. The 
approach allows us to evaluate the excitation energy AE 
and the equal-time one-particle propagator G. We found 
that for a sequence of sudden quenches AE is always 
larger than zero. The thermalization breakdown has a 
dramatic impact on the propagator both at finite times 
and at the steady-state. In particular the steady-state 
G exhibits a power-law behavior \x\~^ with different ex- 
ponents at long and short distances. Remarkable we are 
able to express ft as explicit functionals of the switching 
protocol. We found that at long distances (3 = /?m carries 
informations on the initial state but not on the history of 
the switching protocol, and coincides with the exponent 
of a sudden quench. 21 At short distances j3 = /3 s( j is, in- 
stead, both initial-state and history dependent; it is only 



for continuous protocols that memory is washed out and 
/3 s d equals the exponent of a LL in equilibrium. 

The continuous limit of the recursive procedure leads 
to a simple system of differential equations that can 
be solved numerically. For the excitation energy AE 
we proved that there exist no switching protocols of fi- 
nite duration capable to bring a system from a non- 
interacting ground-state to the ground state of a LL. 
It is only in the adiabatic limit that thermalization oc- 
curs. These results may be relevant to design the optimal 
switching protocol that minimizes the total excitation en- 
ergy or the energy a of given q-excitation. For the prop- 
agator G we clarify how the adiabatic limit is attained 
by studying ramp-switching protocols of increasing dura- 
tion T . At weak coupling we derive an explicit expression 
for G{x,t oo ) that is valid for all x and around the 
crossover distance x w 2vT is dominated by a power-law 
| a; — 2uT| -/3 ( x ) with a x-dependent exponent. The for- 
mula further reproduces with remarkable accuracy the 
different power-laws at long and short distances. Our 
predictions could be experimentally confirmed by mea- 
surements in ultracold fcrmionic atoms loaded in optical 
lattices, recently proposed 35,36 as candidates to realize 
highly controllable and tunable LLs. 
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